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We present a framework for investigating the response of conformally-invariant confined 1+1- 
dimensional systems to a quantum quench. While conformal invariance is generally destroyed in a 
global quantum quench, systems that can be described as or mapped to integrable deformations of a 
CFT may present special instances where a conformal field theory-based analysis could provide useful 
insight into the non-equilibrium dynamics. We investigate this possibility by considering a quench 
analogous to that of the Quantum Newton’s Cradle experiment [JQ and demonstrating qualitative 
agreement between observables derived in the CFT framework and those of the experimental system. 
We propose that this agreement may be a feature of the proximity of the experimental system to 
an integrable deformation of a c = 1 CFT. 


The analytical modeling of the out-of-equilibrium be¬ 
havior of systems subjected to a quantum quench — a 
sudden change in the system’s Hamiltonian parameters 
remains a challenging problem in all but a limited 
number of simple cases. Of particular interest are sys¬ 
tems that do not exhibit simple relaxation to a thermal 
state following a quantum quench; while generically we 
expect systems subjected to sudden changes to eventu¬ 
ally thermalize and reach an equilibrium state as a result 
of interactions, the past decade has seen accummulating 
experimental evidence for and theoretical studies on one¬ 
dimensional systems that do not equilibrate to a simple 
thermal state, but that instead appear to retain memory 
of their initial state (for reviews see [2H1])- Integrability 
has been shown to inhibit thermalization [5], and such 
behavior in experimental systems is often attributed to 
the proximity of these systems to an integrable point. 

The connections between systems that are integrable 
and those that are conformally-invariant have been sub¬ 
ject to ongoing investigation. In particular, certain in¬ 
tegrable field theories can be obtained from massive de¬ 
formations of particular models of conformal field theory 
(CFT) [6], rendering the understanding of thermaliza¬ 
tion within a CFT framework a potentially powerful tool 
for testing some of the ideas arising in the study of the 
connections between conformal invariance and integra¬ 
bility [7]. In principle, if it is known how a particular 
integrable model arises as a perturbation at a conformal 
fixed point, conformal perturbation theory (see, e.g., [5]) 
can be used to compute observables of the integrable 
theory up to arbitrary order. While this approach may 
often become computationally cumbersome beyond the 
lowest orders, it raises the question of whether out-of¬ 
equilibrium analyses of certain CFTs may shed light on 
the post-quench behavior of related integrable models. 

In this Letter we address this question by providing 
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an example of a realistic near-integrable quenched sys¬ 
tem that we show exhibits the behavior characteristic 
of a conformally-invariant system. The system that we 
consider is that of the “Quantum Newton’s Cradle” ex¬ 
periment [TJ, in which an effectively one-dimensional sys¬ 
tem of interacting harmonically-confined bosons was split 
into two oppositely-moving momentum groups; following 
this quench, the system failed to demonstrate any ap¬ 
parent thermalization within experimental time scales. 
While some experimental effects, such as the presence 
of a confining trapping potential in the setup, may in¬ 
troduce weak integrability-breaking effects, the system 
is believed to be well-described by the integrable Lieb- 
Liniger model [9J of delta-interacting bosons, and the 
failure of the experimental system to thermalize has been 
attributed to the integrability of this model. 

To motivate the relation to the CFT picture, we note 
that the non-relativistic Lieb-Liniger model can be ex¬ 
actly mapped to the relativistic sinh-Gordon model in 
an appropriate limit. In particular, under this mapping 
the S-matrix and Lagrangian of the two models coin¬ 
cide mm- The sinh-Gordon model is a massive in¬ 
tegrable deformation of a free scalar field Lagrangian, 
and correlation functions in this model can in principle 
be computed order-by-order in a conformal perturbation 
expansion. 

We proceed as if this system were a c = 1 CFT, 
which is an accurate effective description of the Lieb- 
Liniger model in the limit of either low momenta or hard¬ 
core boson interactions (the latter which map to free 
fermions |12|). This amounts to neglecting higher order 
terms in a perturbative expansion of correlation functions 
of the sinh-Gordon model and hence its non-relativistic 
Lieb-Liniger limit. An important issue when truncating 
such a perturbative expansion is whether higher-order 
perturbative effects, which may not qualitatively change 
the behavior in equilibrium, could have significant effect 
in a non-equilibrium setting on the asymptotic (long¬ 
time) behavior of observables. As we show, this does not 
appear to be the case in a qualitative analysis; we com- 
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ment on this and suggest how a quantitative analysis may 
be performed in order to detect potential deviations. 

In the free boson CFT analyzed here, the harmonic 
confinement of the system implies that (up to an overall 
rescaling) there is a full equivalence up to a phase lag 
of half the system’s size between the position-space en¬ 
ergy density expectation value (given by (T tt )) and the 
momentum-space expectation value. The latter is the 
CFT observable corresponding to the momentum distri¬ 
butions observed in the experiment. Although the ex¬ 
perimental setup was in principle not limited to the low 
momenta or hard-core interactions regimes, we show that 
the experimental momentum distributions and this CFT 
observable qualitatively agree. 

Methods for analyzing quenches in a CFT via a bound¬ 
ary state approach were proposed in EM- In these 
construction^ the system, with Hamiltonian H , is pre¬ 
pared at t = 0 in an initial state |'0o); which is an eigen¬ 
state of a different Hamiltonian Hq. In an analytically- 
continued Euclidean version of the theory this initial 
state can be interpreted as a Euclidean “boundary state” 

| B) that encodes the initial conditions of |^>o)- The 
state of the system at any later time is then given by 
| il>(t)) = e~ lHt |'i/’o)j and correlation functions of observ¬ 
ables can be computed in the boundary state | B), fol¬ 
lowed by an analytic continuation to Lorentzian time. 
Since conformal boundary states are non-normalizable, 
the actual boundary state | B) is taken to be at a certain 
RG distance — the “extrapolation length” To — from the 
conformal boundary state \B) CFT . 

The introduction of this extrapolation length gives the 
boundary state a width of 2r 0 : rather than evolving from 
a single Euclidean boundary state, the time-evolution is 
from the center r = 0 of a slab, or strip, whose top and 
bottom boundaries at r = ±ro correspond to the same 
initial state. Real-time t correlation functions are then 
obtained by evolving observables from r = 0 and analyt¬ 
ically continuing r —> —it. Correlation functions in this 
setup in a 2D boundary CFT (BCFT) pj)] can then be 
computed by making use of the conformal transformation 
that maps the strip to the half-plane. 

Spatial confinement — To make contact with the re¬ 
alistic system we modify this formalism to account for 
systems on a finite interval. We introduce the spatial con¬ 
finement by adding boundaries along the spatial direction 
such that the length of the system is now given by L. The 
resulting boundary state geometry is therefore that of a 
rectangle of length L and height 2r 0 (Fig. |T]) which also 
conformally maps to the half-plane. The transformation 
to the right-half plane [2Dj is given by an elliptic Jacobi 
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FIG. 1: The Schwarz-Christoffel tranformation taking the 
right-half plane to the rectangle with chiral and antichiral 
vertex operators inserted on the vertical midline, r = 0 
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where k is the elliptic modulus, k £ [0,1] and I\ i(...) is 
the complete elliptic integral of the first kind. Its inverse 
is a Schwarz-Christoffel transformation m given by the 
elliptic integral of the first kind 
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and that maps a set of designated points y = 
±i,±i/vT^fc on the imaginary ( x = 0) axis to the 
vertices of a rectangle as shown in Fig. |T] with height 
2ro = 2Ki (\/l — fc) / K\ (k) where tq is the extrapola¬ 
tion length of the previous section. The limit of k —> 0 
corresponds to the infinite-height rectangle (strip) and 
k — > 1 is the limit of zero height. The mapping (HI is 
doubly-periodic (i.e. with one period equal to 2 L) in 
the (real) argumeni]^ for r 0 > 0; this is a feature of the 
open reflective boundary conditions that it imposes, and 
as a result observables in this geometry will display pe¬ 
riodic returns to their initial values, though the period 
may change with the number of insertions. This makes 
this choice of boundary state geometry particularly well- 
suited for modeling harmonically-confined systems. The 
special case of the ground-state expectation value, the 
Casimir energy, has periodicity L as the Jacobi elliptic 
functions only appear squared [20] . 

Vertex operator insertions — We introduce excitations 
in this setup by restricting to c = 1 CFT and consider¬ 
ing local vertex operator insertions. We note that the 
formalism described here carries over with minor modi¬ 
fications to c < 1 minimal models through the addition 


2 The modification by an additive constant here from Ea centers 
the resulting rectangle on the origin of the transformed coordi¬ 
nates. 

3 As a result of the continuation to Lorentzian time all arguments 
considered here are real. 
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of screening charges to all correlation functions. We thus 
consider the free boson action 

s = ^Jd 2 x (ymf , ( 2 ) 

where <f> is a bosonic field, which we take to be com- 
pactified on a circle of radius A, <fi ~ <j> + 2irR , which 
we henceforth set to R = 1. Highest-weight states are 
given by the action of vertex operators on the vacuum 
state, |ra,m) = lim z ,z-yoV nm (z,z) |0), where V nm (z,z) = 
V nm (z) <8> Vnfn(z), and the chiral and antichiral vertex 
operators are given respectively as V nm (z ) =: e ,QWI ^ 2 ) : 
and Vnm(z) =: e lOL ™™<l > ( z ) where m is the winding num¬ 
ber and n is the wave number |22| . 

Their holomorphic and antiholomorphic conformal di¬ 
mensions are given by h nm = (n+ f ) 2 and 

Km = ^ {n - f ) 2 • A Luttinger liquid CFT, 

for instance, is obtained by setting the normalization 
g = K in where K is the Luttinger parameter. The 
bosonic field <j> now represents propagating density fluc¬ 
tuations and is related to the dual variable 6 under the 
T-duality transformation </> <-»• 9 and K o 1/A, n o m. 
For the comparison with the experimental data in the 
following section we will set g = K in subsequent calcu¬ 
lations. 

Split-momentum quench — We implement the split 
momentum quench as a boundary state given by a pair of 
chiral V nm (z) and antichiral Vnm(z) vertex operators of 
the compactified free boson (Fig. [l]) together with Dirich- 
let boundary conditions in conjunction with the elliptic 
Jacobi mapping. The opposite-chirality vertex operators 
act to excite the ground state in analogy to the experi¬ 
mental setup of two excited oppositely-moving clouds of 
bosons. In the CFT analogy, each such cloud is repre¬ 
sented as a peak given by the location of the vertex oper¬ 
ator; in reality the clouds have a certain spread, and we 
later discuss how this spread can be accounted for in the 
CFT analysis. The Dirichlet conditions in conjunction 
with the inherent periodicity of the conformal mapping 
are implemented to mimick the harmonic trapping po¬ 
tential of the experimental setup. 

The Dirichlet condition selects the type of boundary 
states allowed P5H251 . which are given by mm 

||D)} = (4**0"* |(n,0)» D (3) 

nEZ 


where <po is canonically conjugate to the zero mode of the 
free boson and takes values in a circle of radius R = 1. 


The normalization (47 tA') _j (2S] is the g-factor (bound¬ 
ary entropy l|29j for the Dirichlet boundary condition for 
the action pi) with g = K. Unlike in the boundary-less 
case, expectation values of primary operators do not in 
general vanish in a BCFT; in the case of the compact¬ 
ified boson the expectation value can be obtained from 
the boundary states above as 

{Vnm{z,z)) D = ~^=e 1 ^ \z-z\ * o (4) 

V Ao 


where A'o = y/4nK. 

The energy density expectation value (ip o \T tt (t, x) \ ipo) 
at time t for the initial state of the split-momentum 
quench is given, upon analytic continuation t — ► it, by 


^ ((D rec | | (T(w) + T(w)) V nm (W ) Vnm (w')| | -Dree)) 

(5) 

where ||D rec )) is the boundary state state p| follow¬ 
ing the conformal transformation to the rectangle, and 
we have used the decomposition of the energy density 
as the sum of holomorphic and antiholomorphic compo¬ 
nents. Coordinates on the rectangle will be denoted by 
w = x + ir and on the half-plane by z = a + irj. Recall 
that it is the Euclidean time coordinate of the stress ten¬ 
sor, T tt , rather than the time coordinates of the vertex 
operators, that is analytically continued to Lorentzian 
time. The coordinates w' = x' + ir' , where t' = 0, de¬ 
note the location of the vertex operator insertion on the 
rectangle. The equivalence of p]) with the time-evolved 
expectation value of the energy density from the given 
initial state can be understood by noting that the right- 
hand side can be formally expressed as a Euclidean path 
integral with an operator insertion. 

The expectation value pj) can be computed by con¬ 
formally transforming both the vertex operators and 
the stress tensor to the half-plane. The vertex op¬ 
erators transform under the conformal transformation 

as primary fields, V nm (w') = ($p-) V nm (z') and 


V nm (w') = (^= 7 -) V nm (z'), whereas the stress 
tensor acquires the anomalous Casimir term, T(w) = 
(tr 2 3» + Mz’M, where {z;w} = 

{ d 2 z/dw 2 \ ‘ 
l dz/dw J 


(dz/dw) 

is the Schwarzian derivative. In the ab¬ 


sence of spatial boundaries, i.e. the infinite strip limit of 
k —> 0, the Schwarzian derivative term is equal to the con¬ 
stant strip Casimir enery. The Casimir term produced by 
the Jacobi elliptic transformation p) for k > 0 is not a 
constant, and it has a significant qualitative effect on the 
energy distribution. Employing the Ward identity on the 
upper-half plane m 


(T(z)V nm (z')Vnrn(z')) 
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we arrive at the expression for ([5]) 

n 2 

<IW( "'“ ,)> = ^7K (d7IF) 
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where a.h. refers to the antiholomorphic part of the ex¬ 
pression, i.e. 2 z, w —► w, and as a result of the 
Dirichlet boundary condition we have set h = h nm = 
h nm = whs, where as before Kf, = AttK, and made use 
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in computing the chiral-antichiral vertex operator 
correlator. We stress that the coordinates 2 in © must 
be read as functions of the rectangle coordinates w, re¬ 
lated via ([I]). 

Since the transformation Q is from the right-half 
plane, the antiholomorphic coordinates z are rotated 
from the usual upper-half plane ones, i.e. z = —z*. Fi¬ 
nally, the Lorentzian energy expectation value (T tt (f,x)) 
is obtained via a Wick rotation w = x + ir —> x + t and 
w = —x + ir —► —x + t. 

We note that while there appear to be four divergences 
in © for all times t > 0 , in fact two of these divergences 
fall outside of the rectangle boundaries at any given time, 
so that there are effectively only two remaining diver¬ 
gences. These divergences oscillate within the confines 
of the system, coinciding twice within each period of the 
full energy distribution ( T t t(t,x )), which is 2 L as a re¬ 
sult of ©• These divergences are a feature of an analysis 
that - despite the conformal transformation to a finite 
geometry - has been carried out in the thermodynamic 
limit. They are a consequence of the divergence of the 
correlation length in the thermodynamic limit: since the 
system that we consider here is finite of length L, the 
divergences are rounded off owing to the effects of finite 
size scaling [30] (see Appendix [A] for regulation scheme). 

The physical picture that emerges from © is that of 
the non-constant Casimir term, A {z(w), w}+a.h., owing 
to the special type of confinement imposed, competing in 
strength with the two oscilating bumps (regulated diver¬ 
gences) of the terms involving the momentum excitations 
given by the vertex operator insertions. The strength of 

these bumps is given by the vertex operators’ confer- 
2 

mal dimension h = dbj, so that the relative strength of 
these momentum packets to the Casimir term increases 
with lower values of the Luttinger parameter K. Since 
decreasing values of K correspond to increasing values 
of the Lieb-Liniger parameter 7 |31j . their strength in¬ 
creases with 7. 

The elliptic Jacobi transformation with Dirichlet 
boundary conditions thus mimicks the behavior observed 
in the experiment iTj, where the two momentum packets 
repeatedly oscillate within the harmonic trap with a pe¬ 
riodicity such that the two packets coincide twice per 
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period. As a result of the harmonic symmetry of the 
setup, in the case of the non-interacting system assumed 
in the CFT analysis a direct comparison is possible be¬ 
tween position-space and momentum-space distributions: 
at any time t p the momentum-space distribution is seen 
to be equivalent (up to an overall scale factor) to the 
position-space distribution at t x = t p + L/2. In the CFT 
analysis we have assumed that the two packets are highly 
localized; a realistic spread in the momentum may be ac¬ 
counted for by shifting the spatial coordinate away from 
the endpoints of the interval and closer to the middle. 

Fig.[2]shows the regulated (T tt ) plots for two such shifts 
corresponding to a 30% spread and a 50% spread re¬ 
spectively at t x intervals corresponding to integer-period 
t p intervals, and at increasing values of the Luttinger 
parameter. We thus see the characteristic behavior of 
the experimental momentum-space distributions of [T], 
shown for comparison for decreasing initial (input) Lieb- 
Liniger interaction strengths in Fig. [3j The red curves in 
the experimental plots are expanded momentum distri¬ 
butions at single periodic times and are the observable 
most closely expected to correspond to the derived dis¬ 
tributions. 

It would be instructive to investigate whether a more 
quantitative analysis could reveal deviations from the 
CFT description. This can be done by placing the CFT 
parameters in correspondence with the analogous exper¬ 
imental ones, such as putting the Luttinger parameter I\ 
in numerical correspondence with the Lieb-Liniger inter¬ 
action strengths |3lj used in the experiment. The rect¬ 
angle width, L , is representative of the horizontal mo¬ 
mentum scale - but a proper fitting of the experimen¬ 
tal plots necessitates knowledge of the original unsealed 
height of the experimental plots, to be compared with 
(Ttt) after rescaling by the appropriate constant to mo¬ 
mentum space values. The extrapolation length To is 
system-dependent, and the ability to fit this parameter 
could in and of itself provide interesting insight into the 
physics that it models. The spread of momentum in the 
initial peaks must also be taken into account for a proper 
fitting, as demonstrated in Fig. [3] 

Discussion — This intriguing qualitative agreement 
between the experimental distributions of [T] and the cor¬ 
responding distributions for the analogous CFT system 
is not expected for a general system following a quan¬ 
tum quench that injects high energy into the system. It 
may be that the experimental parameters are such that 
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FIG. 2: Derived regulated (see Appendix 0 plots for (Ttt) as a function of position x from the analysis presented here for 
increasing values of Luttinger parameter K at t x = § and subsequent periodic intervals (2 Lm). Up to an overall scale factor 
they correspond to momentum-space distributions at t p = 0 and subsequent periodic intervals. Blue: full fitted distribution; 

dashed red: Casimir energy contribution. Top row: separation of vertex operators is Aa; = 0.5, bottom row: Ax = 0.3. 
Parameters used were L = 1 and extrapolation length to ss 0.278, i.e. k = 0.9999. The charge of the vertex operators was set 

to n = 1. In the plots we have set (j>o = 0 in ([ 6 ]). 
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FIG. 3: Left: Experimental projections jT] of the evolution 
of dephased mommentum distributions without 
thermalization (based on observations after 15 oscillation 
periods) versus actual expanded momentum distributions for 
different dephased peaks coupling strengths 7 d after an 
integer number of periods. From top to bottom: 7 d = 18, 
7 d = 3.2, 7 d = 1.4. The dashed line represents a gaussian 
distribution with the same number of atoms and r.m.s. 
width as the actual distribution. Right: Momentum 
distributions for different coupling strengths obtained by 
transversely integrating absorption images. Top to bottom: 

7 = 4, 7 = 1, 7 =» 0.62. Figures from pQ. 


the c = 1 CFT is still an approximate description of the 
system at the energies used in the experimental setup. 
However, the momenta injected during the quench are 
in principle above those that yield a post-quench Lut¬ 
tinger liquid. The findings of our analysis therefore call 
into question whether special features of the experimen¬ 
tal system — possibly relating to the integrability of the 
system — lead to a post-quench relaxation towards a 
conformal fixed point. 
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Appendix A: Divergence regulation scheme 


In the thermodynamic limit, the correlation length 

_ rp I 

£ ~ f _ ", where t = 1 Tc ° , diverges at the critical tern- 
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peraturc^T = T c . The size of the critical region is then 
given by t ~ ». In a finite-size system the correlation 

length is limited by the system size; the expected scaling 
in a trap of size L is £ ~ L e [32], where 9 is the trap crit¬ 
ical exponent. Experimental systems of trapped ultra¬ 
cold bosons in optical lattices in one dimension are well- 
described by the Bose-Hubbard Hamiltonian [33] and for 
that model it is given by 9 = in the case of a 

power-law potential. This implies that the size of the 
critical region is given by t ~ L~ *. For divergences oc¬ 
curring at x = Xq we therefore place a cutoff at the height 
corresponding to the left boundary of the critical region, 
x = Xq — where a is an arbitrary but consis¬ 

tent choice of constant (a = 0.1 in Fig. [4]), and round 
off the divergences at the corresponding height by find¬ 
ing a best-fit function (skewed exponential ansatz) for a 
set of representative points such that a smooth choice is 
ensured for a given choice of K. We set v = 1 and p = 2 


(harmonic potential) for the distributions derived here. 
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FIG. 4: Fitting scheme for divergence regulation. Left: 

K = 2, Right: K = 4; Ax = 0.5. Curves shown are 
unregulated (T«) plot (solid blue), Casimir hump (dashed 
red), location of singularities (dashed orange), lower bound 
on critical region (dashed purple), corresponding vertical 
cutoff (dashed black straight line), data points used for 
fitting (green), and numerically fitted plots (dashed black) 
using a skew-normal distribution. 
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